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Abstract 

Based on the order parameter expansion, we present an approximate method 
which allows us to reduce large systems of coupled differential equations with diverse 
parameters to three equations: one for the global, mean field, variable and two 
which describe the fluctuations around this mean value. With this tool we analyze 
phase-transitions induced by microscopic disorder in three prototypical models of 
phase-transitions which have been studied previously in the presence of thermal 
noise. We study how macroscopic order is induced or destroyed by time independent 
local disorder and analyze the limits of the approximation by comparing the results 
with the numerical solutions of the self-consistency equation which arises from the 
property of self-averaging. Finally, we carry on a finite-size analysis of the numerical 
results and calculate the corresponding critical exponents. 

Key words: Disorder, Quenched noise. Diversity, Phase-transitions, Global 
coupling. Collective dynamics. Order parameter expansion 



1 Introduction 



The effect of time-dependent noise in extended dynamical systems has been 
the subject of intensive study in the last years [I]. Besides the expected disor- 
dering role, it has been found that some kind of order at the macroscopic level 
can appear by increasing the intensity of the noise. Examples of this paradox- 
ical result include stochastic resonance [2]|3] . or enhancement of the effect of 
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an external forcing under the right amount of noise, coherence resonance jl] 
(also named as stochastic coherence [5]) where a dynamical system displays 
optimal periodicity at the right noise value, noise sustained patterns, struc- 
tures and fronts PI7] . phase transitions where a more ordered phase appears 
when increasing the noise intensity j8|l9] . etc. 



In a very general framework, it has been argued that the resonance with an ex- 
ternal forcing can also be achieved when the time- dependent noise is replaced 
by a more general source of disorder. This includes natural diversity or het- 
erogeneity, competitive interactions, disorder in the network of connectivities, 
etc. and can appear in driven bistable and excitable systems [T0|ll|12] . in lin- 
ear [13] and chaotic [H] oscillators and in a variety of other systems [T5|16|17|18|19|20|21f22f23] ). 
A unifying treatment of the role of noise and diversity for non-forced excitable 
systems, has been developed in 



In this work we examine the effect that structural disorder or diversity, in the 
form of quenched noise, has on some prototypical models of phase transitions 
which have been thoroughly studied in the presence of noise. From the prac- 
tical point of view, the models we will be considering bear some similarities 
with random-field, or impurities, models. As tool of investigation we will refine 
a previously developed order parameter expansion method of approximating 
large systems of coupled differential equations P5|26f27f28f29] with diverse 
parameters. This allows the reduction of the large set of differential equations 
to just three: one for the global, mean field, value and two which describe the 
fluctuations around this mean value. Within this approximation we will ana- 
lyze three different models and show its ease to deliver some understanding of 
the emergent properties of the global behavior. To find the limits of the order 
parameter expansion method we will compare the results with the solution of 
the self-consistency equation which arises from the property of self-averaging. 



The chosen models are a set of globally coupled {/)'*-systems both in the pres- 
ence of additive and multiplicative quenched noise and the canonical model 
for noise- induced phase transition p[9] . It will be seen that quenched noise 
can induce phase transitions (in, out and reentrant) of ordered phases. 



The rest of the paper is organized as follows: In the next section we will de- 
scribe the analytical methods, self-consistency and order parameter expansion; 
in section [3] we will apply those methods to the showcase models and compare 
with the results of numerical simulations; the last section closes the paper with 
the discussion of the presented results. 



2 



2 Models and method 



The type of models we will be considering in this work are defined via differ- 
ential equations for the dynamics of a set of real variables xf. 

Xi = fixi,r]i;X) , i = l,...,N. (1) 

The time derivative Xi{t) = dxi{t)/dt depends on the constant parameter T^j, 
a kind of quenched noise. The set of values {rji, . . . ,rj^} are independently 
drawn from a probability distribution g{ri) of mean H and variance cr^. Cou- 
pling between the different dynamical equations is provided by the presence 
in Eq. ([1]) of the global variable or mean value X{t) = {xi{t)) = j^Y.iLiXiit)- 
For a given realization of the ?7j's variables the Xj's tend in the limit t — )■ oo 
to some asymptotic, stationary values which, in general, will depend on initial 
conditions. Some insight can be obtained if we write Eq. ([T]) as a relaxational 
dynamics [30] in a potential V{xi, rji] X) = — dx/ f{xi', rji] X): 

dV{xi,r]i;X) 
dx, • 

If the potential V{xi, rji] X) is monostable for a particular value of X, then the 
variable Xi{t) tends during the dynamical evolution towards the single mini- 
mum of V{xi, rji] X). Note that the location of this minimum will change with 
time as X evolves. If, on the contrary, V{xi,rii\X) presents several minima, 
the dynamics will tend towards one of the local minimum of the potential. 

In the following we will be interested in characterizing the stationary solution 
by the ensemble average value and fluctuations with respect to realizations of 
the quenched noise and initial conditions of the global variable X. We flrst 
review briefly the self-consistency method and then explain the approximate 
method based on the order parameter expansion. 

2.1 Self- consistency 

This method uses ideas borrowed from the Weiss molecular fleld theory |3T] . 
which is known to be exact for systems with long-range interaction or, equiva- 
lently, in which the interaction occurs through the global variable X, a mean- 
field scenario, as it is our case. Let us denote by x* the stable stationary 
solution of Eq. ([1]). This is nothing but the absolute minimum of the poten- 
tial V{xi,rii\X). It will be a function of rji and the global variable X, i.e. 
X* = x*{T]i, X). For a given realization of the quenched noise variables r^j's, the 
value of the global variable must be obtained from the self-consistency relation 
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X = ■^YliLiX*{r]i, X). It is clear that for large, the sum can be replaced 
by an integral over the distribution g{ri) of the independent T^j's variables: 

X = J dr]g{r])x*{r],X). (3) 



It is then assumed that one can identify the value of X obtained solving this 
equation, as the desired ensemble average, the self-averaging property [32j. In 
general, the possible solutions X of the self-consistency equation ([3]) have to 
be found numerically. A possible scenario is that by changing some parameter 
(e.g. the root-mean square a or the mean H) of the distribution g{T]), the 
solutions bifurcate and the system then presents a phase transition between 
the possible solutions. We will present in the examples below the results of 
this procedure, but will not give any further details about the (in general, very 
involved) numerical method used to solve Eq. ([3]). 



2.2 Order parameter expansion 



For the development of this approximate method we assume, as in the previous 
subsection, that the number of degrees of freedom N is very large and then 
it is possible to substitute the mean value of the distribution g{ri) by the 
system average H = {r]i) = j^Ylf=iVi^ the variance by cr^ = {{r]i — {rii)Y) = 
jfJ2j^i{Vi ~ {VilYi ^iid similar expressions for other cases. 

Our goal is to find an approximate equation describing the dynamics of the 
mean value variable X. To this end, we will expand the evolution equations in 
the deviations ej(t) = Xi{t) — X{t) of the dynamical variables from the mean 
value, and the deviations 5i = rji — H oi the parameters from their mean value. 
The Taylor expansion of Eqs. ([1]) around the mean values up to second order 
gives: 



Xi = /(X, H- X) + e, /,(X, H- X) + 5, /,(X, H- X) + 

^ fUX, H- X) + eA U{X, H;X) + ^ /,,(X, i7; X) + . . . (4) 



, etc. We now take 



With the usual notation /^(X,if;X) = ^^^^^ 

averages and use that (ej) = ^ Yl!i=\ = and = Z^ili = 0. Further- 
more we have {6f) = cr^ as the parameter distribution's variance. So when we 
average over Eq. (jlj) we are left with: 



X = f{X, H; X) + ^f,,{X, H; X) (e^) + /,,(X, H; X) (eA) + 
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f,,{X,H;X) + 0{{e^),{e%), 



(5) 



The evolution of X is then coupled to that of the second moment of the 
snapshot probability density fl = (ef ) = Y^iLi and the so-called shape pa- 
rameter [28j W = {ei6i) = J2f=i ^i^i- We will now obtain evolution equations 
for these two variables. We follow closely the method of [26j but keeping all 
terms up to second order in and 5i. We start by subtracting ([5]) from (jl]) to 
obtain q = — X, which can then be replaced in = (2eiej), W = {5iki). 
After some algebra, and neglecting terms of order 0((ef), {ef6i), . . .) or higher, 
we get: 

X = f{X, H- X) + ^fUX.H- X) + H-X)W + 

2 

^f,,{X,H;X), (6a) 

n = 2n ux, H-x) + 2w h- x), (6b) 

W = W MX, H; X) + a'f,{X, H; X). (6c) 



In summary, within this approximation we have obtained a closed set of three 
differential equations f lBall^ . They have the feature of being coupled only in 
one direction, i.e. W(t) is independent of the others and Q(t) depends only on 



W{t). Steady states with W = Q = X 



OaieW = -a^k- and = a^^'' 



The equilibrium of variable X is given by the solution of: 



fx 



J X 



2 ■ 



J r]ri ~ J XX n2 
J X 



f q 2 f^ "^ f"^ 



fx 



(7) 



where we have simplified notation / = /(X, H] X), etc. As before, an analysis 
of the bifurcations of this equation will allow us to find the possible phase 
transitions of the model. 

Our results, Eqs.dH]), differ slightly from those in the cited sources. In P^27|28] 
the authors require the parameter to be additive, thus setting = 1. In other 
works [26|29] any parameter dependence is allowed, but a coherent regime is 
required, such that terms of order 0{{ef)) and higher are neglected. 



3 Examples 



After presenting the general development of the order parameter expansion 
method, we will now apply it to a few models of relevance in the field of phase 
transitions. Our purpose is to compare the results of our approximation with 
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those coming from the self-consistency equation analysis as well as with numer- 
ical simulations of the different models. Solving the self-consistency equation 
requires in practice a complicated numerical calculation, while our treatment 
is simple and predicts in some cases the existence of phase transitions with 
reasonable accuracy. 



3.1 Globally coupled Landau- Ginzhurg model with additive quenched noise 



The Landau- Ginzburg or 0^ scalar-field has been thoroughly studied from the 
analytical and numerical points of view as a paradigmatic model undergoing 
a second-order phase transition [53]. Here we are interested in this model 
in the case that the stochastic thermal fluctuations have been replaced by 
additive quenched noise, as an example of a random-field scalar model 
The dynamical equations for the set of Xj, i = 1, . . . , A^, real variables are: 

Xi = axi — xl + C {X — Xi) + r]i . (8) 



The study of the model using the self-consistency relation Eq. ^ can be found 
in [12] . Here we want to use the order parameter expansion to derive the main 
properties of this model, in particular the existence of a phase transition as a 
function of the intensity a of the fluctuations of the random fields rji. 

Following the steps from section [21 we obtain the set of equations for the order 
parameter X and the fluctuations VL: 

X = (a-Sn) X - X^ + H (9a) 
tl = 2n{a-C-3X^)+2W (9b) 

W = W [a-C -3X^) +a\ (9c) 
The steady state for the order parameter, Eq. ([7]), leads to: 

0=(a-3 :.^X-X^ + H. (10) 



We now consider the case of zero average field H = (rji) = 0. In that case, 
Eq. ffTOj) can have up to five real solutions. The trivial solution X = 0, always 
exists and it is stable (if C > a) whenever a > ac, with 

rO ifa<0, 
- j - a) if a > 0. ^^^) 
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It turns out that for C > 7a > 0, the set of Eqs. fl9atl9cp contains two additional 
stable fixed point real solutions ±Xq for a < ac- At a = Oc it is Xq = and 
hence Oc identifies a second order, continuous, phase transition (see right panel 
of Fig. [T]). If7a>C>a>0 the range of existence and stability of these two 
additional solutions extends up to cr < ctq, where ctq > Oc is given by: 



^0 = J^(2a + Cf, (12) 



Hence, in the range a G [cXc, o"o] there is bistability between the X = and the 
±Xo solutions. Moreover, two additional symmetric unstable solutions ±Xi 
appear in this range. Therefore, the point ctq signals the appearance of a first 
order, discontinuous, phase transition (see Fig. [1], /e/t). In that range, the three 
stable solutions coexist with the two unstable solutions. 

From a microscopic point of view, the phase transition from the |X| > to 
the X = states can be explained as follows: for a = 0, it is ?7j = 0,Vi; all 
variables end up in the same stationary value = i/a or xi = —y/a and 
the average value satisfies |X| = y/a > . As the noise intensity increases, 
(T > 0, the average value |X| tends to zero and the chances that individual 
values Tji are smaller than —CX grow. This changes the minimum's sign in 
the (individual) potential. As a consequence the distribution of {xi} becomes 
bimodal and the mean value approaches zero. 




0.5 1.0 1.5 2.0 2 4 6 8 10 

a o 



Fig. 1. Bifurcation diagram of the Landau-Ginzburg model with additive quenched 
noise. Order parameter expansion (thick black lines) predict a second order transi- 
tion for C > a {right: C = 10, a = 1) while bistability (first order phase transition) 
appears for C < 7a {left: C = 1.5, a = 1, the unstable solution is plotted as a dotted 
line). The self-consistency solution (grey line) does not show bistability in any case. 
Symbols show the results of numerical simulations of the evolution equations aver- 
aged over 10^ realizations of the quenched noise variables r]i and initial conditions. 
N = 10^,10^,10^ (circles, squares, diamonds, respectively). 

The existence of a phase transition from order to disorder predicted by the 
order parameter expansion simple approximation scheme is confirmed by the 
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numerical solution of the self-consistency equation ([3]) [12]. However, the tran- 
sition appears to be always second-order, so indicating the validity of the pre- 
diction of the approximate order parameter expansion in the limit of large 
coupling. In fact, the critical value cTc predicted by the order-disorder transi- 
tion, Eq. f lTTj) . deviates systematically from that coming from the numerical 
integration of the self-consistency equation ([3]) for large values of the coupling 
constant C, as shown in figure [2l although the relative error between the two 
values decreases as C increases. 




C 



Fig. 2. Critical intensity of the additive quenched noise for the Landau-Ginzburg 
model versus coupling strength for a = 1. Prediction of order parameter expan- 
sion (llip as continuous line, exact solution ^ as dotted line. The order parameter 
expansion predicts a bistability region for C < 7a. 

We have also compared these predictions versus the results coming from inten- 
sive numerical simulations. In the simulations we have integrated the full set 
of equations ([H]) up to the steady state and, then, we have computed the order 
parameter m = {{\X\)) and its fluctuations X = ^ [((^^)) ~ (d-^l))^]- Here 
X = Y^iLi and ((■ ■ ■)) denotes an ensemble average with respect to realiza- 
tions of the random variables rji and initial conditions. The simulation results 
for the order parameter are indicated by symbols in figure [TJ As usual, the 
transition from order to disorder is smeared out due to finite-size-effects but 
the numerical simulations do approach the results of the self-consistency equa- 
tion as the number of variables increases. We have analyzed our data using 
standard finite-size-scaling relations [35)136] and found that the dependence 
of the order parameter on a can be well fitted by m{a, N) = N'^f^ fmi^N^) 
with e = 1 — o"/(Tc, 6 ~ 0.33 and being a scaling function, see evidence in 
the left panels of figure ([3]) for two different values of the coupling constant. 
Note that this scaling relation implies that in the thermodynamic limit, the 
order parameter vanishes as m{a) ~ {a^ — o")^^^, the typical mean-field result. 
Similarly, the fiuctuations can be fitted by the form x((T, A^) = N'^f^ {eN^^ , 
with c ~ 0.67 and the appropriate scaling function, as demonstrated in the 
right panels of figure ([3]) again for two different values of the coupling con- 
stant. This implies that in the thermodynamic limit, the fluctuations diverge 
as ~ Ic^c — o"!^''' with 7 = c/h ~ 2. 
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Fig. 3. Finite-size scaling analysis of the Landau-Ginzburg model with additive 
quenched noise. Rescaled simulation data for low coupling (top graphs, Oc = 1.094) 
and high coupling {bottom graphs, = 6.203). Ensemble average m {left) and 
fluctuations xi'^wht) as defined in the main text. Exponents: b = 0.33, c = 0.67. 
Ensemble sizes: N = 10^10^,10^ (circles, squares, diamonds). In all cases: a = 1. 
Here, Cc has been determined to a high degree of accuracy by using the numerical 
solution of Eq. ([3]). 

3.2 Globally coupled Landau-Ginzburg model with multiplicative quenched 
noise 



We now consider the case in which the quenched noise couples multiplicatively 
to the variable xf. 

Xi = {a + Tji) Xi — xl + C {X — Xi) . (13) 



This model has been extensively studied in the case that the T^j's are indepen- 
dent white noises and it has been found that an increase in the noise intensity 
leads to a transition from disorder to order p!|37f38f39] . We want to compare 
the predictions of the self-consistency equation with the order parameter ex- 
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pansion and numerical simulations to check if a similar result holds in the 
case of quenched noise. Without coupling (C = 0) Eq. ( !T3|) is a prototype of 
supercritical pitchfork bifurcations (see e.g. in |1D]) with two possible sets of 
solutions: = is the stable solution whenever a + ?7j < or Xj = ±-y/a + r]i 
are stable solutions and = is unstable for a + rji > 0. 

To study the consequences of coupling, C > 0, we use the above developed 
order parameter expansion approximation. After setting H = (rji) = 0, the 
equations are: 

X = aX-X^- 3Xn + W, (14a) 
tl= (2a- QX^ -2C)n + 2XW, (14b) 
= (a - - C) + _ ^^^^^ 

The equilibrium condition ([7]) leads to: 



3X^2 Xa^ 
{C-a + 3X^f ^ C-a + 3X2 



= aX - X=^ - — , 2 + 7^ ^ , ■ (15) 



Similarly to the uncoupled case this equation has two different regimes of 
solutions: On one hand, if a > the stable solutions of Eq. f|T^ are X = ±^/a 
for 0" = 0. As (T increases, |X| monotonically increases as well (see fig. HI left). 
On the other hand, if a < then X = is a stable solution for small a. At 
some value Uc it becomes unstable and a fork of solutions grows out of zero 
(see fig. m right). Oc is determined by Eq. flT^ and is related to a and C by: 

al = a{a-C). (16) 



Of. identifies a second-order phase transition from disorder to order (i.e. from 
X = to X 7^ 0). In this case of a < 0, the value cXc grows monotonously with 
coupling strength C, a rather counterintuitive observation, since it means that 
the coupling hinders the ordering and more structural disorder is needed to 
induce macroscopic order (fig. |5]). 

The numerical solution of the self-consistency equation ^ is qualitatively 
similar to the results of the order parameter expansion approximation, how- 
ever |X| doesn't increase monotonically with increasing a. It rather reaches 
a maximum and decreases after that approaching zero asymptotically. Note 
that his is not a (reentrant) phase transition since |X| = is only reached for 
(J — )► oo. 

The simulation results for the order parameter are shown as symbols in figure 
m At this scale no finite-size-effects can be seen at the phase transition. In a 
thorough data analysis with finite-size-scaling relations at cTc, in the way we 
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Fig. 4. Bifurcation diagram of the Landau-Ginzburg model with multiplicative 
quenched noise. Positive values of a show order without noise {left: a = 0.1), 
whereas negative value show order only with a finite value of the noise intensity 
{right: a = —0.5). In both panels it is C = 10. The order parameter expansion ap- 
proximation scheme gives a monotonous solution while the exact solution of Eq. ([3]) 
reaches a maximum and decreases for large a (grey line). Symbols are the result 
of direct numerical simulations of Eqs. ()13p averaged over 10^ realizations of the 
quenched noise variables rji for A'^ = 10^, 10^, 10^ (circles, squares, diamonds). 




I 

5 10 15 20 

c 

Fig. 5. Critical noise for bifurcation versus coupling strength for a = —0.5 for the 
Landau-Ginzburg model with multiplicative quenched noise. The prediction by the 
order parameter expansion is shown as continuous line, the exact solution ([3]) is 
shown as dotted line. 

did in the first example, we found exponents of 6 ~ c ~ 0.5 to fit the order 
parameter and fluctuations (see fig. E]). These scaling relations imply, again in 
the thermodynamic limit, that the order parameter vanish and the fluctuations 
diverge as m(a) ~ (cTc — 0")^^^ and x(cr) ~ |(Tc — cr\^'^, with 7 = c/b ^ 1 
respectively. 



3.3 Canonical model for noise-induced phase transitions 



At the last example we will study a model for which a genuine phase transition 
induced by multiphcative noise has been shown pi9] with the feature that the 
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Fig. 6. Finite-size scaling analysis of the Landau-Ginzburg model with multiplicative 
noise for a = —0.5, C = 10. Rescaled ensemble average m (left) and fluctuations 
X (right) of 10'^ numerical simulations with = 10^,10^,10^ (circles, squares, 
diamonds). Critical point, as from Eq. (I3|), is Uc = 2.169; exponents: b = 0.5, c = 0.5. 

ordered phase is reentrant, it only exists for intermediate noise intensities. The 
equation for an individual element is: 



[l + xfj" + [l + xfj r,, + C{X 



Xi J 



;i7) 



and the reduced system according to section I2l2] (again setting (r/j) = 0) reads: 

,2 1 



X = -X 1+X 



12 1 + X^ X-8X 



n + 2XW (18a) 



-{l + X'^y -AX^{l + X'^^+C +2W{l + X 
{l+X' 



n = 2n 
w = w 

The equilibrium condition ([7]) becomes 



AX^ 1+X^] +C 



(7^ ( 1 + 



(18b) 
(18c) 



-X (l + X' 



2 (- (6 + 6X2) X- 4X3)^2 



2n2 



+ 



;i + 6X2 + 5X4 + C)' 

Xa2(l + X2) 



+ 



+2 



1 + 6X2 + 5X4 + C 



(19) 



Equation ( 1191) has the stable solution X = for a < cTc or a pair of symmetric 
solutions X 7^ for 0" > ctc. Here X = becomes unstable (see fig. [7j). The 
value of (Tc indicates the location of a second-order phase transition. It follows 
from analyzing the Jacobian of f ll8aH18cj) and calculates to: 



a. 



l + C 
V2C-4 



(20) 
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Accordingly, a minimal coupling C > 2 is necessary to induce the phase tran- 
sition. An analysis of this relation shows that cJc has a minimum with respect 
to C. Therefore, see figure [HI the transition is predicted to be reentrant with 
respect to C: the ordered phase only exists in a range of values for C, with the 
surprising prediction that too a large coupling destroys the ordered state. The 
predictions of the order parameter expansion are in qualitative agreement with 
those obtained after solving the self-consistency equation. However, whereas 
the order parameter expansion predicts incorrectly that the order parame- 
ter monotonously increases with a, as shown in figure [TJ the self-consistency 
equation instead predicts that the transition is reentrant also with respect to 
the quenched noise intensity cr, see figure [HI Both reentrant behaviors were 
observed in the case of time-dependent noise [8|9j . 

a 




Fig. 7. Bifurcation diagram of model (|17p {left). Order parameter expansion (thick 
black line) and exact solution (grey line) together with the ensemble average of 10^ 
numerical simulations for N = 10^,10^,10^ (circles, squares, diamonds). Coupling 
is C = 10. On the right the unsealed fluctuations are shown. 




Fig. 8. Critical noise intensity versus coupling strength for model (|17p. The predic- 
tion ()20p of the order parameter expansion (continuous line) and the exact solution 
([3]) as dotted line. Only the latter shows reentrance with respect to a. 



Again we have compared the predictions with the numerical integration of 
the set of equations (fTTj) . The simulation results are shown as symbols in 
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figure [71 Due to finite-size-effects the theoretical results are approached with 
increasing number of particles, reentrance and the dependence of cTc from 
C are observed. Analyzing the data as we have done with the other examples, 
we find exponents for the scaling relations of 6 ~ 0.33 and c ~ 0.67. As in the 
first case this implies the relations m{a) ~ (ctc — cr)^^^ and x{^) ~ — o"| 
7 = c/6~2in the thermodynamic limit. Figure |9] summarizes the fitted 
simulation data. 
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Fig. 9. Rescaled simulation data for the first {top) and the second (bottom) phase 
transition (C = 10). Mean value (left) and fluctuations (right) for ensembles of 
N = 10^,10^,10^ (circles, squares, diamonds). Critical points, as from Eq. are 
(7c = 2.749 for the first and Cc = 5.169 for the second transition. Exponents are: 
6 = 0.33,c = 0.67. 



4 Conclusions 



In this paper we have constructed an approximate analytical scheme based 
on the order parameter expansion [25|26f27f28f29] to study the macroscopic 
behavior of extended systems which are globally coupled. We have used the 
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method to study in detail the phase diagram of three widely used models of 
phase transitions in scalar systems: the Landau-Ginzburg scalar model with 
both additive and multiplicative quenched noise and a genuine model for noise- 
induced phase transitions where time-dependent noise has been replaced by 
quenched, time-independent noise coupled multiplicatively to the dynamical 
variable [8]. 

We have compared the results of our simple approach with those coming from 
a numerically involved, but in principle exact, treatment based on the self- 
consistency relation and with extensive numerical simulations of the corre- 
sponding dynamical equations for each model. In the case of additive noise, 
the main result is that macroscopic order is destroyed when increasing the 
intensity of the quenched noise. In the other two cases, when noise appears 
multiplicatively, we find that macroscopic order appears for intermediate value 
of the intensity of the quenched noise. Since the quenched noise can represent, 
for instance, diversity or heterogeneity, it appears paradoxically that some 
amount of structural disorder is needed in order to observe macroscopic order. 

Furthermore it has been shown numerically, that all investigated models follow 
a finite-size scaling law and the exponents have been determined. It suggests 
a common universality class for the Landau-Ginzburg model with additive 
quenched noise and the canonical model for noise induced phase transitions, 
whereas the Landau-Ginzburg model with multiplicative quenched noise yields 
different exponents. A more detailed analysis of the finite-size relations and 
their possible dependence with the system parameters will be presented else- 
where |41j . 

The method of order parameter expansion, which we lead consistently up to 
terms of second order, is a tool which reduces large systems to only a couple of 
reduced variables. The advantage is its very easy management. In this paper 
we have proven that reliable conclusions can be drawn with that method in 
some cases. Since it is an expansion around mean values the method yields 
good results for low values of the intensity of the quenched noise or for high 
synchronization of the subunits. Otherwise, the method might not be reliable. 
As a consequence the reentrant phase transitions were not predicted in the 
studied cases for multiplicative noise. It is an open issue how to modify the 
method in order to predict the reentrant transitions. 



5 Acknowledgements 

The authors acknowledge financial support from the EU NoE BioSim, LSHB- 
CT-2004-005137, and project FIS2007-60327 from MEG (Spain) and FEDER 
(EU). NK is supported by a grant from the Govern Balear. 



15 



References 

[I] J. Garcia-Ojalvo and J.M. Sancho. Noise in Spatially Extended Systems. 
Springer-Verlag, New York, 1999. 

[2] L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni. Stochastic resonance. 
Rev. Mod. Phys., 70:223, 1998. 

[3] P. Hanggi and F. Marchesoni, eds. Topical issue on stochastic resonance. In 
The European Physical Journal, volume 69, 2009. 

[4] A. Pikovsky and J. Kurths. Coherence resonance in a noise-driven excitable 
system. Phys. Rev. Lett, 78:775, 1997. 

[5] A.A. Zaikin, J. Garci'a-Ojalvo, R. Bascones, E. Ullner, and J. Kurths. Doubly 

stochastic coherence via noise-induced symmetry in bistable neural models. 
Phys. Rev. Lett, 90:030601, 2003. 

[6] M. Santagiustina, P. Colet, M. San Miguel, and D. Walgraef. Noise-sustained 
convective structures in nonlinear optics. Phys. Rev. Lett., 79(19):3633-3636, 
Nov 1997. 

[7] M.G. Clerc, C. Falcon, and E. Tirapegui. Front propagation sustained by 
additive noise. Phys. Rev. E, 74(1):011303, Jul 2006. 

[8] C. van den Broeck, J.M.R. Parrondo, and R. Toral. Noise-induced 
nonequilibrium phase transition. Phys. Rev. Lett., 73:3395, 1994. 

[9] C. van den Broeck, J.M.R. Parrondo, R. Toral, and R. Kawai. Nonequilibrium 
phase transitions induced by multiplicative noise. Phys. Rev. E, 55:4084, 1997. 

[10] C.J. Tessone, C.R. Mirasso, R. Toral, and J.D. Gunton. Diversity-induced 
resonance. Phys. Rev. Lett, 97:194101, 2006. 

[II] C.J. Tessone, A. Scire, R. Toral, and P. Colet. Theory of collective firing induced 
by noise or diversity in excitable media. Physical Review E, 75, 2007. 

[12] R. Toral, C.J. Tessone, and J.V. Lopes. Collective effects induced by diversity 
in extended systems. Eur. Phys. J. Special Topics, 143:59-67, 2007. 

[13] R. Toral, E. Hernandez-Garcfa, and J.D. Gunton. Diversity-induced resonance 
in a system of globally coupled linear oscillators. International Journal of 
Bifurcations and Chaos, 19:3499, 2009. 

[14] H. Chen and J. Zhang. Diversity-induced coherence resonance in spatially 
extended chaotic systems. Physical Review E, 77, 2008. 

[15] M. Gosak. Cellular diversity promotes intercellular ca2-|- wave propagation. 
Biophysical Chemistry, 139:53, 2009. 

[16] E. Ullner, J. Buceta, A. Diez-Noguera, and J. Garcia-Ojalvo. Noise-induced 
coherence in multicellular circadian clocks. Biophysical Journal, 96:3573, 2009. 



16 



[17] D. Zanette. Interplay of noise and coupling in heterogeneous ensembles of phase 
oscillators. European Physical Journal B, 69:269, 2009. 

[18] C.J. Tessone and R. Toral. Diversity-induced resonance in a model for opinion 
formation. European Physical Journal B, 71:549, 2009. 

[19] S. Postnova, K. Voigt, and H.A. Braun. A mathematical model of homeostatic 
regulation of sleep- wake cycles by hypocretin/orexin. Journal of Biological 
Rhythms, 24:523, 2009. 

[20] H.S. Chen, Y. Shen, and Z.H. Hou. Resonant response of forced complex 
networks: The role of topological disorder. Chaos, 19:033122, 2009. 

[21] D. Wu, S.Q. Zhu, and X.Q. Luo. Cooperative effects of random time delays 
and small-world topologies on diversity-induced resonance. European Physics 
Letters, 86:50002, 2009. 

[22] M. Perc, M. Gosak, and S. Kralj. Stochastic resonance in soft matter systems: 
combined effects of static and dynamic disorder. Soft Matter, 4:1861, 2008. 

[23] J. A. Acebron, S. Lozano, and A. Arenas. Amplified signal response in scale-free 
networks by collaborative signaling. Physical Review Letters, 99:128701, 2007. 

[24] C.J. Tessone, A. Scire, R. Toral, and P. Colet. Diversity-induced resonance. 
Phys. Rev. E, 75:016203, 2007. 

[25] S. de Monte and F. d'Ovidio. Dynamics of order parameters for globally coupled 
oscillators. Europhys. Lett, 58:21-27, 2002. 

[26] S. de Monte, F. d'Ovidio, and E. Mosekilde. Coherent regimes of globally 
coupled dynamical systems. Phys. Rev. Lett, 90(5):054102, 2003. 

[27] S. de Monte, F. d'Ovidio, Hugues Chate, and E. Mosekilde. Noise-induced 
macroscopic bifurcations in globally coupled chaotic units. Phys. Rev. Lett., 
92:254101, 2004. 

[28] S. de Monte, F. d'Ovidio, Hugues Chate, and E. Mosekilde. Effects of 
microscopic disorder on the collective dynamics of globally coupled maps. 
Physica D, 205:25-40, 2005. 

[29] I. Gomes Da Silva, S. de Monte, F. d'Ovidio, R. Toral, and C. R. Mirasso. 
Coherent regimes of mutually coupled chuas circuits. Physical Review E, 
73(036203), 2006. 

[30] M. San Miguel and R. Toral. Stochastic effects in physical systems. 
In J. Martinez, E. Tirapegui and R. Tiemann, editors. Instabilities and 
nonequilibrium structures VI, pages 35-120. Kluwer academic publishers, 2000. 

[31] H.E. Stanley. Introduction to phase transitions and critical phenomena. Oxford 
university press, 1971. 

[32] D.P. Landau and K. Binder. A guide to Monte Carlo simulations in statistical 
physics. Cambridge university press, 2000. 



17 



[33] D.J. Amit and V.M. Mayor. Field Theory, the Renormalization Group and 
Critical Phenomena. World Scientific Publishing Co.Pte. Ltd., 3rd edition, 
2005. 

[34] A. P. Young, editor. Spin Glasses and Random Fields. World Scientific 
Publishing Co.Pte. Ltd., 1998. 

[35] J.L. Cardy. Finite-Size Scaling. Elsevier science publishers, 1988. 

[36] H.P. Deutsch. Optimized analysis of the critical behavior in polymer mixtures 
from monte carlo simulations. J. Stat. Phys, 67:1039, 1992. 

[37] C. Van den Broeck, J.M.R. Parrondo, J. Armero, and A. Hernandez-Machado. 
Mean field model for spatially extended systems in the presence of multiplicative 
noise. Phys. Rev. E, 49(4):2639-2643, Apr 1994. 

[38] J. Garci'a-Ojalvo, J.M.R. Parrondo, J.M. Sancho, and C. Van den Broeck. 
Reentrant transition induced by multiplicative noise in the time-dependent 
ginzburg-landau model. Phys. Rev. E, 54(6):6918-6921, Dec 1996. 

[39] J.M.R. Parrondo, J. Buceta and F.J. de la Rubia. Random ginzburg-landau 
model revisited: Reentrant phase transition. Physical Review E, 63:031103, 
2001. 

[40] S.H. Strogatz. Nonlinear dynamics and chaos: with applications to physics, 
biology, chemistry and engineering. Addison- Wesley, 1994. 

[41] N. Komin and R. Toral. work in progress. 



18 



